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Abstract 



Viscoelastic and thermodynamic properties of transient gels formed by telechelic polymers are studied on the basis 
of the transient network theory that takes account of the correlation among polymer chains via network junctions. The 
global information of the gel is incorporated into the theory by introducing the elastically effective chains according 
to the criterion by Scanlan and Case. We also consider effects of superbridges whose backbone is formed by several 
chains connected in series with several breakable junctions inside. Near the critical concentration for the sol/gel 
transition, superbridges becomes infinitely long along the backbone, thereby leading to the short relaxation time r of 
the network. It is shown that r is proportional to the concentration deviation A near the gelation point. The plateau 
modulus Goo increases as the cube of A near the gelation point as a result of the mean-field treatment, and hence 
■^j- | the zero-shear viscosity increases as ijq ~ GooT ~ A 4 . The dynamic shear moduli are well described in terms of the 
Maxwell model, and it is shown that the present model can explain the concentration dependence of the dynamic 
moduli for aqueous solutions of telechelic poly (ethylene oxide). 
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I. INTRODUCTION 

Transient gels formed by associating polymers have been attracted widespread interests in recent years. [l| As- 
sociating polymers are polymer chains carrying specific groups capable of forming aggregates through noncovalent 
bonding. 0, H 0, H, H, 0, H, H EH EH Above a certain polymer concentration, they form a transient gel by connecting 
sticky groups on polymers. This transformation is thermoreversible in general. In the previous paper[l2j (referred to 
as I in the following), we presented a theoretical model of transient gels formed by junctions comprised of limited num- 
ber of hydrophobic groups with an intention to understand thermodynamic properties of linear rheology of telechelic 
associating polymer systems. As the first attempt, elastically effective chains (or active chains) were defined locally, 
i.e., chains whose both ends are connected to other end groups are elastically effective irrespective of whether these 
groups are incorporated into an infinite network (gel) or not. We could explain, to some extent, the concentration 
dependence of the dynamic shear moduli described in terms of the Maxwell model, but the sol/gel transition of the 
system could not be treated properly due to the local definition of active chains. 

In this paper, we take account of the global information of the infinite network by making use of the criterion 
suggested by Scanlan |l3| and Case [l4[ for a chain to be active. This criterion states that telechelic chains are 
elastically effective if their both ends are connected to junctions with at least three paths to the infinite network. We 
assume that these chains deform according to the macroscopic deformations applied to the gel. Static properties of 
transient gels have been studied by Tanaka and Ishida on the basis of this criterion. [l6| We here consider not only 
primary active chains (referred to as primary bridges in this paper) but also active superchains (called superbridges) 
whose backbone is an aggregate of several bridges connected in series. Effects of superbridges cannot be negligible 
especially when one study dynamical properties of transient gels because they enhance the relaxation time of the 
network due to a number of internal junctions possible to break as suggested by Annable et aZ..[2| We can describe 
the transition between the sol state and gel state (appearance of the infinite network) in this theoretical framework. 
The critical behavior of viscoelastic quantities near the sol/gel transition point is shown to be much affected by 
superbridges. 

It has been established up to now that telechelic polymers self-assemble in dilute solution to form flowerlike micelles. 
Pham et al. @, 0] indicated that the solution of flowerlike micelles resembles a colloidal dispersion of adhesive hard 
spheres in the concentration dependence of the shear modulus. Quite recently, Meng and Russel showed that 
the colloidal theory describing the nonequilibrium structure of dispersions under shear explains the high-frequency 
plateau modulus of telechelic poly(ethylene oxide) (PEO). In this paper, we attempt to theoretically describe linear 
rheology of telechelic polymers in the absence of intramolecular associations. It can be shown that experimentally 
observed dynamic shear moduli that are characterized by the high-frequency plateau modulus, zero-shear viscosity 
and relaxation time, are well described in terms of this theoretical treatment. It indicates that the transient network 
theory is useful tool not only for the study of rheological but also for thermodynamic properties of transient gels 
when it is extended so that both the correlation among polymers and global structure of the network are taken into 
consideration. 

This paper is organized as follows. In section [Til we will review assumptions and definitions employed in I. They 
are also employed in this paper. In section HTT1 linear viscoelasticities of the transient gel will be studied within the 
framework of the Scanlan-Case criterion for active chains. As the first step, only primary bridges are taken into 
consideration in this section. Effect of superbridges will be discussed in section IIV1 In section |Vj liner rheology 
including effects of superbridges will be studied. Section IVT1 will be devoted to a summary. 

II. ASSUMPTIONS AND DEFINITIONS OF FUNDAMENTAL QUANTITIES 

We consider athermal solutions of n linear polymers (or primary chains) per unit volume. Functional groups capable 
of forming junctions through noncovalent bonding are locally embedded in both ends of the primary chain. Common 
assumptions employed in this paper and in I are as follows: 1) any number of functional groups are allowed to be 
bound together to form one junction; 2) reactions among functional groups are allowed to occur only in a stepwise 
manner; 3) primary chains are Gaussian chains; 4) the Rouse relaxation time of the primary chain is much smaller 
than the characteristic time of the macroscopic deformation applied to the system and the lifetime of association 
among functional groups; 5) the looped chain formed by a single primary chain is absent; 6) the molecular weight 
of the primary chain is much smaller than the entanglement molecular weight, and hence effects of the topological 
interaction among chains is ignored. 

The terminology used in this paper (and I) are as follows: 1) the number of functional groups forming a junction 
(or the aggregation number), say fc, is referred to as the junction multiplicity; 2) the junction with the multiplicity 
k is called the fc-junction; 3) the primary chain whose head is incorporated into the fc-junction whereas whose tail is 
a member of the fc'-junction is referred to as the (fc, fc')-chain (we virtually mark one end of each chain to identify a 
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head and tail of the chain for convenience, although physical properties of the chain is homogenous along the chain); 
4) the primary chain whose one end, irrespective of whether it is a head or tail, is incorporated into the fc-junction is 
called the fc-chain. 

As in I, we define F ktk i (r, t)dr as the number of (k, fc')-chains at time t per unit volume with the head-to-tail vector 
r ~ r + dr. The total number of (k, fc')-chains (per unit volume) is then given by 



V k .k'{t)= JdTF k , k ,{T,t)=V k ,, k (t) (II.l) 

where two subscripts of v k ,k' {t) are exchangeable due to the symmetry along the polymer chain. The total number of 
primary chains 

« = EE ( n - 2 ) 

fc>l k'>l 

is conserved independent of time. The number of fc-chains is given by 

Xk(t)=Y,^,k'(t), (H.3) 

k>>l 

and then we can express the number of fc-junctions as 

/**(*) = CM 

The number of functional groups belonging to fc-junctions is kji k (t) = 2xk{t) while the total number of functional 
groups is 2n, so that the probability that an arbitrary chosen functional group is in a fc-junction can be expressed as 

q k (t) = (H.5) 
n 

Eq. (|II.2|) is equivalent to the normalization condition of q k , i.e., 

E«*(*) = 1 - ( IL6 ) 

fe>i 

The extent of reaction (or the probability for a functional group to be associated with other groups) is written as 

aW-E*w = 1 -^w- ( IL7 ) 

k>2 

III. THEORY ON THE BASIS OF THE SCANLAN-CASE CRITERION FOR ACTIVE CHAINS 

A. Global Structure of the Network 

Here, we briefly review how to treat the global structure of the network according to refs. 

EE E3, Ell- Above 

a certain concentration of primary chains, an infinite network (gel) is formed. In the postgel regime, the extent of 
reaction a' with regard to the functional groups in the sol part is different from that a" in the gel part EMEU- Eq. 
(III. 71) is the average extent of reaction for all functional groups in the system, i.e., it can be expressed as 

a(t) = a'(t)w s (t) + a"{t)w G {t) (III.l) 

where ws(t) is the fraction of the sol part and wa(t) = 1 — Ws(t) is that of the gel. The sol fraction can be written 

as El El E3 

w s (t)=J2lk(t)^, (HI.2) 

k>l 



where Co is the probability that a randomly chosen unreacted group belongs to the sol part through its main chain. 
In the pregel regime, we have only Co = 1- In the postgel regime, on the other hand, we have Co less than 1. Thus Co 
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elastically effective chain 

path to the gel network 
dangling end 



/ 

O'=4,fc=6)-junction (;'=3,£'=4)-junction 

FIG. 1: The classification of the junction by the multiplicity k and the path connectivity i to the gel network. For example, 
a left-hand junction is formed by six functional groups but it is connected with the gel network only through four paths 
(represented by arrows); the other two paths are connected with dangling ends (indicated by dotted circles). A primary chain 
whose both ends are connected to the junctions with i > 3 is elastically effective. 




is useful as an indicator of gelation. A primary chain belongs to the sol if its both ends are associated with the sol 
part, so that the sol fraction can be also expressed as 

«sW = (E«*c*)Co _1 ) a - (IIL3) 
fe>i 



Therefore, Co is the root of the following equation 
where 



x = u(x), (III.4) 



u(x) = J2Qk(t)x k -\ (III.5) 



k>\ 

If (|III.4|) has more than one root, we must employ the smallest one. 

Now we consider the connectivity of a functional group to the gel network according to the theoretical treatment 
by Pearson and Graessley [l7|. Let fi^f. be the number of junctions with multiplicity k that is connected to the gel 
network through i paths (0 < i < k). Such a j unction is called the (i, fc)-junction in the following. According to the 
multinomial theorem, it takes the form [l6|, [lTj 

u\ 

= MfcW .^.j CrT'a - Co) 1 - (m.6) 

Then the number of paths emersed from (i, fc)-junction is Xi,k{t) — (i/2)fXi i k(t). Here we employ the criterion of 
Scanlan [l3[ and Case [3] to decide whether the primary chain is active or not. They suggested that the primary 
chain whose both ends are connected to junctions with the path connectivity larger than or equal to 3 is elastically 
effective (see FigHJ) . According to this criterion, the number of active chains whose one end is belonging to fc-junctions 
is 

k 

xf(t) = £>,*(*) = X k(t)(l - Co)[l - Co" 1 - (fc- l)C fe - 2 (l - Co)]- (IIL7) 

Note that x\ (*) = X% f (*) = 0- The total number of active chains is then given by [l6| 

^ e// W = £%f (*) = "(1 - Co) 2 (l - u'(Co)), (III.8) 

fc>3 

where we have used a relation Co = u (Co)- The number of active (fc, /c')-chains (i.e., the chain whose one end is 
connected to the (i, fc)-junction with i > 3 and the other end is belonging to the (i', fc')-junction with i' > 3) can be 
defined as 

The following relation 



£ uf k ,(t)=^(t) (IIL10) 



k,k'>3 



holds as it should be. 
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B. Time-Development of Chains 



The number of (k, fc')-chains with the head-to-tail vector r obeys the following time-evolution equation: 
dF k<k ,(r,t) 



dt 



+ V-(r fc)fc ,(r,*)F M ,(r,i)) = W k)k >(r,t) (for k,k' > 3), (III.ll) 



where T kik i(r,t) is the rate of deformation of the head-to-tail vector r of the (k, &/)-chain. When a macroscopic 
deformation is applied to the gel, only active chains deform. Some (k, fc')-chains are active but the other (k, fc')-chains 
are not because each junction has the different path connectivity even if it has the same multiplicity. To take this 
into account, we put 

ik,k'(t) =Pk,k>(t)k(t)r (111.12) 
where k(t) is the rate of deformation tensor applied to the gel, and 

p ktk ,(t) = j^y- IH.13 

is the probability for a (fc, fc')-chain to be active. Eq. (|III.12[) states that active chains deform affinely to the 
macroscopic deformation on average. Eq. holds only for k,k' > 3 because both 1-junction and 2-junction 

cannot have path connectivity more than or equal to 3 and chains connected with such junctions do not deform. We 
assume that these elastically ineffective primary chains are Gaussian chains; i.e., the probability distribution function 
that they take the head-to-tail vector r is 

fo(r)=(^X /2 e Xp (-i^) (111.14) 



\2%Na 2 J ^\ 2Na 2 J 

where N and a is the number and the length of the repeat unit constructing a primary chain, respectively. For 
example, the distribution function for the dangling chain is written as F kt i(r,t) = Vk,i(t)fo(r). The right-hand side 
of pil.lip represents the net increment of the (fc, fc')-chain per unit time due to the reaction between the end groups 
on the (k, fc')-chain and the groups on the other chains. As shown in I, it takes the form 

W Kk , (r, t)=- [/3 fc (r) + [3 k> (r) + B k {t) + By (t) + P k (t) + P v (t)] F k . v (r, t) 
+ [p k -i(t)v ltk , (t) + p fe ,_ x (i)i/ M (t)] / (r) 
+B k+1 {t)F k+1M . (r, t) + B k . +1 (t)F fcl fc/+i (r, t) 

+P k -i(t)F k -i, k '(r,t) +P k >-i(t)F k>k >-i(r,t) (for k,k' > 3) (111.15) 

where f3 k (r) is the probability that a functional group is dissociated from the fc-junction per unit time, and p k (t) is 
the probability for an unreacted group to be connected with the fc-junction per unit time. As in I, we assume that 
the dissociation rate does not depend on the junction multiplicity nor the cnd-to-cnd length of the primary chain that 
is connected to the junction (P k (r) — [3). On the other hand, the connection rate p k of the functional group with a 
fc-junction is assumed to be proportional to the volume fraction of the fc-junction, and we put p k {t) = f3\ipq k (t)h k 
where A = exp(e/ksT) is the association constant (e is the binding energy for the attraction of functional groups), 
tp = 2nvQ is the volume fraction of the functional group (i>o is the effective volume of a segment), and h k is the factor 
depending on k that give a limitation on the junction multiplicity. In the following, we use c = \ip = 2\<j)/N as a 
reduced polymer concentration (0 = Nnvo is the volume fraction of the primary chain). Under these assumptions, 
B k (t) and P k (t) in (|III.15|) are expressed as B k = (3{k — 1) and P k (t) = f3kcqi(t)h k , respectively. 
The kinetic equation for (k, fc')-chains has the same form as the one derived in I, i.e., 

= w ktk ,(t)+wk>,k(t) (for k,k' > 1), (111.16) 



(11 



where 



Wk,k>(t) = -/3fc(l + cqi(t)hk)vk,k'(t) + 0kUk+i,k'(t) + (k- l)(3cqi(t)hk-ivi<.-i,k>(t) 

+pch k -.iq h - 1 (t)v 1 ,k>(t) (forfe>2), (III.17a) 

wi,v{t) = p(j2 p w(t) + V2,k>(t)) - (3c(j2 h m{t) + h iqi {t)y iM (t). (III.17b) 



l>2 l>l 
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Note that (jlll.lip reduces to (|III.16|) by integrating it with respect to r (for k,k' > 3). The kinetic equation for 
fc-chains is obtained by taking a summation in 1 6[) for overall k'. It takes the form 

= v k (t) (for k > 1), (III. 18) 

at 

where 

v k (t) =-/3k(q k (t) - q k +i(t)) + /3cfc(/i fc _ig fe _i(t) - h k q k {t)) qi (t) (for k > 2), (III.19a) 

+ <&(*)) - 0c(j2hm(t) + h iqi (t)) qi (t). (III.19b) 

l>2 1>1 

Once q k (t) is derived by solving (|III.18j) . we can obtain the number of (fc, /c')-chains from the relation Vk. k '{t) = 
nq k (t)q k '(t) and ( from pil.4[) . respectively. Then v kk ,(t) and P k , k '(t) = ^f{,/(£)/^fc,fc'(*) can be obtained. By 
putting (t) into (jlllllj) and solving the equations, we can derive F kik > . 

We study two special models of junctions by putting a limitation on the multiplicity, i.e., the saturating junction 
model and the fixed multiplicity model [H, EH, EH- In the saturating junction model, we allow junctions to take only 
a limited range of multiplicity k = 1,2, ...s m . This condition is realized by putting h k = 1 for 1 < k < s m — 1 and 
h k = for k > s m . If s m = oo, junction can take any value of multiplicity without limitation. In the fixed multiplicity 
model, all junctions can take only the same multiplicity s (except for k = 1). This situation can be approximately 
attained by introducing a small quantity 5(<C 1) and by putting h k = 5 for 1 < k < s— 1, = <5 _ ' s_2 ' ) , and h k = 
for k > s — 1.1121 We set S = 0.01 in most cases. 



C. Equilibrium Properties 

By putting dq k /dt = in (|III.18[) . we can obtain q k in equilibrium state. It turns out to be 

q k - JkC^ql (111.20) 

where ^ k = Yi^i hi f° r k > 2 and ji = 1. The fraction of unreacted groups q\ is determined from the normalization 
condition 7(2)91 = 1, where j(z) = J2 k >i r JkZ k ~ 1 and z = cq%. The function defined by (|III.5j) can be expressed as 
u(x) = r y(zx)/ r y(z) in the equilibrium state, and hence £0 is a solution of the following equation for a given z (or c): 

» = 2M. (111.21) 

7(2) 

In the case that the junction can take any value of multiplicity without limitation, for example, z is smaller than 1 
(see I) and hence 7(2) = 1/(1 — z) (we can put h k = 1 for all k in this case). Thus (0 is obtained as the smallest root 
of the equation x = (1 — z)/(l — xz), i.e., 

Hv-i (°- s <^)' <IIL22) 

where z* = 1/2. Note that z* is interpreted as the critical value of the parameter z for the sol/gel transition. The 
number of elastically effective chain given by (|III.8|) is then analytically expressed as a function of z: 

eft ( (0 < z < z*) , 

Vn = {n(2 { z-z*),zf (2*~<2<1) ' ( IIL23 ) 

We can also express these quantities as a function of c by the use of the relation z — c/(l + c), that is, 

Co = {l /c %t*) ' (IIL24) 



and 



' eff = \ ( ( *w^! C<C *! ; (HI.25) 
1 n((c - c*)/c) (c> c*) y 
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(a) (b) 




FIG. 2: (a) Reduced sol/gel transition concentration for the saturating junction model (circles) plotted against the maximum 
multiplicity s m and for fixed multiplicity model (triangles) plotted against the multiplicity s. (b) Sol/gel transition curves 
drawn in the temperature-volume fraction plane for saturating junction model with Ao = 1 and TV = 100. The maximum 
multiplicity is varying from curve to curve. See also Fig|5] 



where c* = 1 is the critical reduced concentration for gelation (see Figf5] (a) for the case of limited multiplicity). 
Near the sol /gel transition point, the number of elastically effective chains increases as the cube of the concentration 
deviation. [16j| i.e., v e ^ ~ A 3 where A = (c — c*)/c*. In the high concentration limit, on the other hand, all primary 
chains become elastically effective, i.e., v e ^ — > n for c — > oo. 

In general, the sol/gel transition point is obtained as the point at which Co becomes lower than 1. This is equivalent 
to the point where the weight-average molecular weight of the cluster diverges. For polycondensation by multiple 
reaction, Fukui and Yamabe have shown that this condition (an appearance of the macroscopic cluster) is given 
by@ 

(f w - 1)(m™ - 1) = 1 (HI.26) 

where f w {=2 for the present model) is the weight-average functionality of the primary chain, and /i w is the weight- 
average multiplicity of the junctions given by 

fi w = V^ fe = l + ^M. (111.27) 

1{Z) 

A boundary curve separating the sol region and gel region in the temperature-concentration plane can be drawn 
by the use of the equation 2\(T*)<p* /N = c* where T* and (f>* are the critical temperature and volume fraction 
of the primary chain, respectively, and c* is obtained according to the procedure stated above. Since the binding 
free energy is comprised of the energy part eo and the entropy part So, the association constant is rewritten as 
A(T) = Aoexp(Xo/T) with Ao = exp(— Sq/Rb) and To = ea/ks- Thus we have 

T* = T / log (Nc*/2\q4>*) . (111.28) 

FigEl (b) shows the sol/gel boundary lines given by pil.28[) for the saturating junction model as an example. 



D. Dynamic-Mechanical and Viscoelastic Properties 



We now apply a small oscillatory shear deformation with an amplitude e to the gel network. A xy component of 
the rate of this deformation tensor is represented by K xy (t) = eu cos cot while the other components are (lo is the 
frequency of the oscillation). Let us expand F^ t y (r, t) with respect to the powers of e up to the first order: 

F kM , (r, t) = Fi% (r) + eF$, (r, t). (111.29) 

We can put (see I) 

F k%( r ) = "k,k'fo(r), (III.30a) 

= (ffU.Msin^ + 5 ^,Hcos^)^/o(r). (III.30b) 
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The number of (k, /c')-chains does not depend on time as far as the small shear deformation is concerned, [20, [2T], [22| 
so that Vk,k'(t) and Qk(t) can be represented by their equilibrium values Vk,k' an d <7fc given above, respectively. Then 
the probability for a (k, fc')-chain to be elastically effective is given by its equilibrium value Pk,k' — y kk'l Vk k' wu "h 

eff _ [1 - Ct 1 ~ (fc ~ l)Cp fc - 2 (l - Co)] [1 - Co*'' 1 - (k> l)Co fc '- 2 (l - Co)] mj „ n 

" k ' k ' ~ ^ l- Z1 '{C z)h(z) ■ (UL31) 

The in-phase g' kk , (lu) and out-of-phase g kk >(Lu) amplitude of F kk , are directly related with the storage and loss 
modulus of (k, fc')-chains through the relations G' kk ,(uj) — kBTg k k ,(ui) and G' kk ,(Lu) = k B Tg k k ,(u)), respectively. 
Since primary chains are symmetric along the chain, two subscripts are exchangeable, i.e., g' k k , — g k , k (and g kk , = 
9k' k)- Note that g' k 1 = g k 1 = g k 2 = g k2 = for k > 1 since 2-chains and 1-chains are effectively in equilibrium 
state. By substituting (|III.29[) with pil.30[) into , we have a set of equations for g' k k , and g k k , as follows 

9k,k> = (-Qk,k'9k,k> + B k+ i9' k +i,k' + Bjf+iSfcV+i + Pk-i9k-i,k> + Pk'-i9'ik'-i) I" + < f ,l, (III.32a) 
9k,k> = (Qk,k'9kM' - Bk+ig'k+i.k- - B k'+i9k,k'+i - p k-i9k-i,k' - p k'-i9k,k'-i) lu (III.32b) 

(for k, k' > 3) 

where B k = f3(k - l),Pk = (3kcqih k , Qk,k'(t) = (3k(l + cq\hk) + (3k' {I + cqihy), and u*jff k , is given by (|III.31[) . It 
should be emphasized here that the last term in the right-hand side of (|III.32ap is v kk , instead of Vk,k' ( see I as a 
reference). This is a consequence of the assumption (|III.12[) . 

The total moduli within the framework of the Scanlan-Case criterion for elastically effective chains are given by 
taking a summation of G' k k , (u>) and G' k k , (lu) over k, k' > 3, i.e., 

G'H - k B T £ £ g' k>k , [to), (III.33a) 

fc>3 fe'>3 

G"(u) = k B Tj2J29k,k'(u). (IIL33b) 

fc>3 fc'>3 

They are well described in terms of the Maxwell model with a single relaxation time (not shown here). In the high 
frequency limit, (|III.32a[) becomes g' k k , (lu — > 00) = v kk , . Therefore, the plateau modulus defined by = G'(lo — > 00) 
can be expressed as 

G x = V 4fk B T = nk B T{\ - Co) 2 ( 1 - Z J^f] . (111.34) 

The reduced plateau modulus GaoKuksT) coincides with the number of active chains derived by Tanaka and Ishida 
for telechelic polymers. As they have shown, it agrees well with experimental data for aqueous solution of hy- 
drophobically modified ethylene oxide- urethane copolymers (called HEUR) reported by Amiable et aL.[2| On the 
other hand, the relaxation time r of the gel determined from the peak position of pil.33b[) does not agree well with 
experiments as in the case of I because it depends on c only weakly. This discrepancy can be ascribed to the absence 
of superbridges in elastically effective chains. A superbridge is a linear cluster of primary chains whose backbone 
includes a number of junctions with the path connectivity i = 2 to the gel network. Both ends of a superbridge are 
connected to junctions with the path connectivity i > 3. Fig[3] shows an example of the superbridge whose backbone 
is formed by four primary chains. In the Scanlan-Case criterion, superbridges are not regarded as the elastically 
effective chains; only primary bridges, or active primary chains (see FigJ3]), are assumed to be responsible for the 
elasticity of the network. If we take superbridges into account, not only that the number of active chains becomes 
larger, but that the relaxation time of the network becomes shorter since their lifetime is shorter than that of primary 
bridge due to breakable node inside the backbone. We will discuss the effects of superbridges on G'(u>) and G"(u>) in 
the next section. a 



a We use the nomenclature superbridge after ref. where a cluster formed by connecting flower micelles linearly through bridges is called 
the superbridge. Their effects on the relaxation time is roughly estimated in this reference. In this paper, we count the number of 
superbridges in more detailed way by making use of the path connectivity to study their effect on the relaxation time. 



FIG. 3: Upper: schematic representation of a superbridge (surrounded by dotted ellipsoid). Lower: schematic representation 
of a primary bridge (or active primary chain) . 



end primary chain 




internal primary chain 



FIG. 4: A schematic of a superbridge (only paths to the network are drawn). It is comprised of two end primary chains and 
several internal primary chains. The number of end primary chains is written as m(i > 3, i' = 2) whereas the number of internal 
primary chains is m(i = 2, i' = 2). 



IV. EFFECTS OF SUPERB RIDGES 



A. Number of Superbridges 

Let m(i, i') be the number of primary chains whose one end is connected to the junction with the path connectivity 
i to the gel whereas the other end is connected to the junction with the connectivity i'. The number of primary 
bridges is then represented as b 

m(i > 3,i' > 3) = n(l - Co) 2 (l - £ ^p) = 4c- ( IV -!) 
The total number of primary chains incorporated into the gel through both ends is 

m(i >2,i'> 2) = ]T]T;tt,fc = n{\ - Co) 2 = i< eff , (IV.2) 

k>2 i>2 

and the number of primary chains comprising backbones of superbridges is 

m(i >2^ = 2) = = n(l - Co) 2 ^ 1 ^ = < eW - (IV-3) 

fc>2 1{Z) 

A relation v e ^ — Vg C + v^ eud holds as it should be. Near the sol/gel transition concentration (or A= (c — c*)/c* <C 1), 

these quantities increase as v^ eud ~ v e ^ ~ A 2 (and <~ A 3 ) because 1 — Co is proportional to A while zj'(( z)/j(z) 
is proportional to c. Since the superbridge is comprised of two end primary chains and a number of internal primary 



b We denote the number of primary bridge as v^q instead of u e ^ in the rest of this article. 
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FIG. 5: The total number of elastically effective chain (solid line), the number of primary bridges or elastically effective chains 
defined on the basis of the Scanlan-Case criterion (dotted line), the number of super bridge (broken line), the number of primary 
chains per a super bridge (dash-dotted line), the ratio between the number of superbridges and the number of primary bridges 
(dash-double dotted line) plotted against relative concentration deviation A = (c — c*)/c* for the saturating junction model 
with the maximum multiplicity fixed at s m — 15. 



chains (see Fig0]), we have a relation v e ^ seud — m(i = 2, %' = 2) + m(i > 3, i' = 2) where 

(eff \2 

/• o •' o\ -eff X2,k X2,k' \ U pseud) nxr a\ 

m(i = 2,i = 2) = ^ ^ v "—jf — r = (IV.4) 



k>2 k'>2 

is the number of internal primary chains, and 



eff eff 



m(i > 3, i> = 2) = vf seud m(i = 2,i' = 2)= % e/ f e " rf (IV.5) 

is the number of end primary chains. In (|IV.4|) . X2,k/v e ^ is the probability for a chain in the network to be connected 
with the (i = 2, £;)-junction, and hence v e ^ (X2,k/^ )(X2,k' /v e *f) is the number of primary chains whose one end is 
connected with the (i — 2, fc)-junction while the other end is belonging to the (i' = 2, fe')-j unction. The number of 
superbridges is a half of the number of end primary chains of superbridges, i.e., 

eff eff 

v su P er = I ,. > 3 / = 2 ) = VSC ' V P f Seud , (IV.6) 

Therefore, the total number of elastically effective chain is turned out to be 

, .eff _ ,eff , super 

--^-^O-^^w)- (,v - 7) 

FigE] shows vfj tal together with the number v^q of primary bridges, the number i/ su P er of superbridges, and the 
relative amount of the superbridges v su P er Jv e g C compared with primary bridges as a function of A for the saturating 
junction model (s m = 15). The number of superbridges increases {v su P er A 3 ) near the sol/gel transition concentra- 
tion, but it decreases at higher concentration since the number of dangling ends decreases. Thus a peak appears in 
v super a ^ a jnodegt concentration. It should be emphasized that v super increases with decreasing A, and it finally 
reaches 0.5 for A — > although both v su P er anc l become close to in this limit. It indicates that the effect of 
superbridges cannot be ignored as compared with primary bridges especially in the vicinity of the sol/gel transition 
point. Fig also show the number of primary chains forming a superbridge given by 

_ jpse^^ (Iy g) 

ySuper v ' 

With decreasing A, many primary chains become incorporated into superbridges as Sb ~ 1/A (for A « 1) indi- 
cating that the superbridge becomes longer along the backbone. With increasing concentration, on the contrary, Sb 
approaches to 2 since m(i = 2, i' = 2) becomes close to 0. Summarizing, in the vicinity of the sol/gel transition 
concentration, (i) the superbridge is comparable in number to the primary bridge although both are few, and (ii) the 
superbridge is infinitely long. 
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B. Breakage Rate of Superbridge 

Let us here focus on a primary chain whose one end is connected to the junction (say A) with the path connectivity 
*A > 3 while the other end is belonging to the junction (A') with the path connectivity ix> > 2. Such a primary 
chain is elastically effective. If ix> > 3, then the chain is a primary bridge, so that the dissociation rate of the end 
group from the junction A' is (3. If «a' = 2, on the other hand, the chain is the end primary chain of the superbridge 
(see FiglH]), and we assume that the breakage rate of internal chains from the junctions B,B',C,... is reflected in the 
dissociation rate of this end primary chain from the junction A'. Then we can put it as the sum of its own dissociation 
rate /3 and the dissociation rate 2(s(, — l)/3 of 2(sb — 1) functional groups B,B',.-- on the internal primary chains. As 
a result, the dissociation rate of A' on average can be expressed as (3 + 2(s;, — l)p(3 = (3 e " , where p is the probability 
for i A . to be 2 and is given by p = m(i >3,i' = 2)/(m(i > 3, i' > 3) + m{i >3,i' = 2)) = vf seu J(v e ^ + vf seud ). We 

replace (3 in (|III.32p with (3 e ^ in the following in order to incorporate the short lifetime of superbridges into account. 
It should be noted that the equilibrium condition pil.20[) still holds after this replacement, so that the discussion 
given in IIII CI and IIV Al of this paper and IV in I remains valid. When A is small, (3 e ^ is inversely proportional to A 
since Sb ~ 1/A while p ~ 1 for A <C 1. Therefore, we see that the relaxation time r of the gel that is approximately 
given as the reciprocal of (3 e " is proportional to A near the sol/gel transition concentration. 

Let Pffl? 1 be the probability for a (k, fc')-chain to be a primary bridge or an end primary chain of a superbridge. 
It is given by 

/ eff \ 

P, total \ total/h,!*' fyxr n\ 

k,k> — ~ ) l iV - 9 J 



Vk,k' 



where 



, eff ( eff -eff -eff eff \ 

(vf t ,)k k> = ^77^77 + ^77^T (IV.10) 

V total/ *,* 2 \ ,,eff yeff jjeff eff J V ) 

\ v sc v sc J 

is the number of (k, fc')-chains that is the primary bridge or the end primary chain of the superbridge. (Xk = 
Si=2 = Xfc(l — Co)(l — Co _1 ) i s tne number of paths (> 2) emersed from the k(> 2)-junction, and hence xf 
is the probability that a fc-junction satisfies the condition i > 2. On the other hand, X^k l v sc ' s ^ ne probability for 
a fc-junction to fulfill the condition i > 3. Eq. pV.lOp is expressed in a symmetric form.) Eq. pV.10[) satisfies the 
following relation: 



EE = pv.ii) 

fc>2 fc'>2 

as it should be. d If a macroscopic deformation is applied to the gel, not only primary bridges but also superbridges 
deform accordingly. According to pil,12[) . we assume that the rate of deformation vector of the (fc, fc')-chain is given 



, e // \. .. _ ,tff 



c Internal primary chains of the superbridge (i.e., chains whose both ends are connected to the junctions with the path connectivity 
two) are also elastically effective in a sense that they are ingredients of the superbridge. However, these chains can be treated in 
this theoretical framework only indirectly as in the previous theory. 0| Effects of these internal chains are taken into account via the 
end primary chains of the superbridge in the rest of this paper. For example, the breakage rate of internal chains is reflected in the 
dissociation rate of the end primary chains. Furthermore, the elasticity of superbridges is represented by that of the end primary chains 
by assuming that the end primary chains deform affinely (see (IIV. 121 0, This is valid when the macroscopic deformation applied to the 
gel is small as in the present case. 

d Note that xf f /"sc = °- 
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FIG. 7: (i) The dynamic shear moduli (reduced by nksT) for the saturating junction model as a function of the frequency. 
The relative concentration deviation A = (c— c*)/c* from the sol/gel transition concentration c* is varying from curve to curve, 
while the maximum multiplicity of the junction is fixed at s m = 20. (ii) The reduced plateau modulus, (iii) reduced zero- 
shear viscosity, and (iv) relaxation time plotted against the relative concentration deviation for several maximum multiplicities 
(increasing from bottom to top). 



by 

v k ,k'(t)=Pl^ l k{t)v. (IV.12) 

By substituting (|IV.12j) into (jlll.llj) . we can obtain the equations for g'S"2 that is given by l|III.32|) but with {vfj ta i)k,k' 

instead of v^l,. Therefore, in the high frequency limit, g' k k , reduces to {v^ tal )k,k' ■ The total (observable) modulus is 
given by 

C<">M = *flT££42( w )' (IV.13) 

k>2 k'>2 

and then we find G' — > G^ — kgTu^^ for ui — » oo. 

V. RESULTS AND DISCUSSIONS 

A. Saturating Junction Model 

FigEI (i) shows the dynamic shear moduli calculated from pV.13[) for the saturating junction model. Relative 
concentration deviation is varying from curve to curve. They are well described in terms of the Maxwell model with 
a single relaxation time. Near the sol/gel transition concentration (A <C 1), the plateau modulus and the relaxation 
time increase as Goo ~ A 3 (see FigJTJ (ii)) and r ~ A (Fig[7J (iv)), respectively. As a result, the zero shear viscosity 
increases as rjo ~ GoqT ~ A 4 (FigJT] (iii)). Note that these powers stem from the mean-field treatment. 6 For example, 
we can explain r ~ A as follows. Let ^ be the radius of gyration of the superbridge. If we assume that the superbridge 



c Rubinstein and Semenov have found the same power laws from the mean-field treatment for multifunctional polymers that can connect 
with each other through pairwise association between functional groups on polymers. |24H 
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FIG. 8: (i) The dynamic shear moduli (reduced by nksT) obtained for the saturating junction model as a function of the 
frequency. The maximum multiplicity s m is varying from curve to curve with a relative concentration deviation fixed at A = 1. 
(ii) The reduced plateau modulus, (iii) reduced zero-shear viscosity, and (iv) relaxation time plotted against the maximum 
multiplicity for several relative concentration deviation (increasing from bottom to top). 



1/2 

obeys the Gaussian statistics, it is estimated to be £ ~ s b . On the other hand, £ corresponds to the network mesh 
size, and hence it obeys the scaling law £ ~ A _l/ with v — 1/2 for A <C 1. By comparing two expressions, we can find 
si, ~ 1/A. Thus the mean lifetime of bridges (primary bridges and superbridges) that corresponds to the relaxation 
time of the network is approximately estimated to be r ~ l/[j3 + 2(s& — l)p/3] ~ A//3. FigJ5]shows the dynamic shear 
moduli as a function of the maximum multiplicity s m . The relaxation time increases with s rn because the number of 
superbridges decreases as s m increases. 

The reduced plateau modulus explicitly depends only on the reduced polymer concentration and the maximum 
multiplicity of the junction. Therefore, it is written as 

-^=fi(c,s m ). (V.l) 

Similarly, the reduced zero-shear viscosity and the relaxation time can be expressed as 

ftno 

nksT 

and 



/ 2 (c,s m ) (V.2) 



/3r = / 3 (c,s m ), (V.3) 

respectively {f\ ~ fa are dimensionless functions of c and s m ). In order to investigate how the dynamic shear moduli 
depend on the temperature and the polymer volume fraction, let us rewrite (|V.1[) ~ (|V.3[) as 

iSi = ^/i( C (T/T o ,iV,0), Sm ), (V.4) 

^ = ^e^Mc(T/T ,N,^s m ), (V.5) 

Pot = e T °/ T -7 3 (c(T/T , N, cf>),s m ), (V.6) 

respectively. We have put W — em the derivation of (|V.5|) and (|V.6j) . In this case, the dissociation rate at temperature 
T is written as 

[3 = co e-^ k - T = poe 1 -^'/? (V.7) 
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FIG. 9: (i) The plateau modulus, (ii) zero-shear viscosity, and (iii) relaxation time plotted against the polymer volume fraction 
for several temperature with s m = 20, iV = 100 and Ao = l. Marked points A, B, C, D on the horizontal axis of (i) indicate the 
critical volume fraction (f>* for each temperature corresponding to the points in Fig[2](b). (iv) Arrhenius plots of the zero-shear 
viscosity and the relaxation time with = 0.05, s m = 20, iV = 100 and Ao = l. 



where /3q is the dissociation rate at temperature To. Recall that the reduced concentration depends on the temperature, 
molecular weight, and the polymer volume fraction as c(T/To, N, 4>) = 20Aoe T °/ T /N . Thus, for example, the zero- 
shear viscosity (|V.5|) depends on the temperature through j3, c and a prefactor. The unitless plateau modulus (|V.5|) . 
zero-shear viscosity (|V.5[) and relaxation time (|V.6|) are shown in Fig[5] (i) ~ (iii) as a function of the polymer volume 
fraction for several temperature. Volume fraction is varying across the sol/gel transition line drawn in Figj2] (b) for 
each temperature. If the volume fraction (f> is small, Goo and 770 depend on <fi through the reduced concentration c 
and a prefactor, while they are approximately proportional to <f> if 4> (and hence c) is large because /1 and fa depend 
only weakly on c (see FigtZJ) in this case. On the other hand, r depends on tf> only through c. As shown in Figj9] (iv), 
the zero-shear viscosity and the relaxation time approximately show the Arrhenius law temperature dependences. At 
higher temperature, we can see a slight deviation from the Arrhenius law. This deviation stems from the fact that r/o 
and t depend on T not only through (3 (see (IV.7|) ) but also through c (and a prefactor in the case of tjq), On the other 
hand, 7/0 and r depend only weakly on c at lower temperature, and hence they show approximately the Arrhenius 
law. We can guess from (|V.5|) and (|V.6p that the dynamic shear moduli at temperature T can be superimposed to 
the curve at the reference temperature T re f if they are horizontally and vertically shifted by a factor of ay and bx, 
respectively, where 



ax = exp 



1 



br 



T ref 

T re f fl{c{T re f/T ,N, 



f 3 (c(T/T o ,N,0),s m ) 
h{c{T re }/To, N, 4>), s m ) ' 

1 1 Sm. ) 



T fMT/T ,N,<l>),a m ) ' 
Especially for larger tj> or lower T, (|V.8[) is approximately written as 



(V.8a) 
(V.8b) 



br 



exp 
T 



-T 



1 



T, 



ref 



(V.9a) 
(V.9b) 



It has been revealed by Annable et al. that the shift factor given by (|V.9[) produces the master curve successfully. @ 
In Fig llOl (i). we compare theoretically predicted dynamic shear moduli (i.e., plateau modulus, zero-shear viscosity 
and relaxation time) with experimental observation for aqueous solutions of telechelic PEO of 20 kg/mol [ll| and 35 
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FIG. 10: Comparison between the theoretically obtained plateau modulus (top), zero-shear viscosity (middle), relaxation time 
(bottom) for the saturating junction model (lines) and experimental data obtained for (i) tclechelic PEO with narrow molecular 
weight distribution and fully end-capped with Cie alkanes (M lII =20,000 for HDU-20 pj and M W =35M0 for HDU-35 Q) and 
(ii) HEUR end-capped with the same alkanes (M„=20,000 for hd-20 and M ro =33,100 for hd-35 0). Values of molecular 
parameters s m , P and A«o used to draw theoretical curves are listed in TABLE H] 

kg/mol [tJ with narrow- molecular- weight distribution and fully end-capped with Ci6 alkanes. We call these polymers 
HDU-20(35) according to ref..[ll| The reduced concentration c was converted into the polymer concentration in weight 
percentage c w through a relation c = (2Q00Na/M)Xvqc w (Na is Avogadro's number). We have three molecular 
parameters for a given molecular weight: s m , \vq and f3. (Note that (3 is not required to calculate Goo.) Values of 
these parameters used to draw theoretical curves are listed in TABLE fl] We find better agreements between theory 
and experiment than in the case that the short lifetime of superbridges is not taken into consideration, [l^] The value 
of A^o increases with increasing the molecular weight. This indicates that the effective volume vq of a functional 
group increases with increasing the chain length. In Fig llOl (ii). we attempt to fit theoretical curves to experimental 
data reported by Annable et al. [2J for HEUR of the similar molecular weight (but with broader molecular weight 
distribution) and end-capped with Cie alkanes. They are called hd-20(35) after ref..[lll Parameter values adopted to 
fit experimental data are also listed in TABLE |TJ We still find a good agreement between theory and experiment in 
spite of a broader molecular weight distribution of hd polymers. A difference in the value of \vq between HDU and 
hd for each (averaged) molecular weight might stem from a difference in the polydispersity of the PEO backbone. A 
ratio between the value of Xvq for HDU-20 and for HDU-35 (3.2) is close to that between hd-20 and hd-35 (2.9). A 
discrepancy in the value of (3 for HDU and for hd might stem from the difference in the coupling agents between the 
alkanes and the PEO backbone as suggested in refs..[7l. [ll| 
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polymer 


$ rn 


P [1/b] 


A« x 10 23 [m 3 ] 


HDU-20 


20 


60 


0.1 


HDU-35 


20 


60 


0.32 


hd-20 


20 


2.3 


0.08 


hd-35 


20 


2.3 


0.23 



TABLE I: Values of molecular parameters used in Fig llOl 
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FIG. 11: (i) The dynamic shear moduli (reduced by nfesT) obtained for the fixed multiplicity model as a function of the 
frequency. The relative concentration deviation is varying from curve to curve, while the multiplicity of the junction is fixed at 
s — 5. (ii) The reduced plateau modulus, (iii) reduced zero-shear viscosity, and (iv) relaxation time plotted against the relative 
concentration deviation for several multiplicities. 

B. Fixed Multiplicity Model 

Fig llll shows the dynamic shear moduli of the fixed multiplicity model together with the plateau modulus, zero- 
shear viscosity and relaxation time plotted against the relative concentration deviation. These quantities obey the 
same critical behavior as in the saturating junction model, i.e., they increase as Goo ~ A 3 , 770 ~ A 4 and r ~ A with 
increasing A near the sol/gel transition concentration. 

In FigIT21 theoretical curves are compared with experimental data for telechelic PEO. Values of parameters used 
to draw theoretical curves are listed in TABLE [TTJ We find a disagreement between theory and experiment in the 
relaxation time; theoretical curves increase more rapidly with increasing the concentration than experimental data, 
and they become approximately flat above a certain concentration. This is because the fraction of junctions with the 
path connectivity i = 2 (or, in other words, the number of superbridges) is small when the junction can take only 
a single multiplicity/ This tendency becomes more pronounced with increasing the multiplicity because the fraction 
of i = 2 junctions is smaller for larger multiplicity. Thus we guess that the multiplicity should not be fixed at a 
single value. In real systems, junctions might be cores of flower-like micelles and the aggregation number (i.e., the 
number of chains per junction), say Sfi ower , is almost independent of the polymer concentration as some researchers 



The relaxation time does not depend on the concentration for the fixed multiplicity model if effects of superbridges are not taken into 
account. 1 12H 
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FIG. 12: Comparison between the theoretically predicted plateau modulus (top), zero-shear viscosity (middle), relaxation time 
(bottom) for the fixed multiplicity model (lines) and experimental data observed for (i) telechelic PEO with narrow molecular 
weight distribution fully end-capped with Cie alkanes (M w =20,000 for HDU-20 [Tl[ and M„=35,000 for HDU-35 @]) and (ii) 
HEUR end-capped with the same alkanes (M„=20,000 for hd-20 and M„=33,100 for hd-35 @]). Values of molecular parameters 
s, P and \vq used to draw theoretical curves are listed in TABLE HT1 



polymer 


s 


a [i/s] 


\v x 10 23 [m 3 ] 


HDU-20 


6 


90 


0.09 


HDU-35 


6 


90 


0.32 


hd-20 


6 


4 


0.07 


hd-35 


6 


4 


0.23 



TABLE II: Parameter values used in Fig[T2l 



indicated. In this case, the number of bridge and dangling chains emersed from the junction (i.e., multiplicity of the 
present theory) can be less than Sfi ower . This situation corresponds to the saturating junction model, not the fixed 
multiplicity model. This might be the reason why the saturating junction model can describe the dynamic shear 
moduli of telechelic PEO better than the fixed multiplicity model. 
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VI. SUMMARY 



We developed the theory of transient networks with junctions of limited multiplicity. The global information was 
incorporated into the theory by introducing the elastically effective chains (active chains) according to the criterion by 
Scanlan and Case and by considering the effect of supcrbridges whose backbone is formed by several chains connected 
in series. Linear viscoelasticities of the network were studied as functions of thermodynamic quantities. Near the 
critical concentration for the sol/gel transition, supcrbridges are infinitely long along the backbone and their number 
is comparable with that of primary bridges. Thus the mean lifetime of bridges is quite short near the critical point 
and so does the relaxation time. It was found that the relaxation time is proportional to the concentration deviation 
A near the sol/gel transition concentration. Since the plateau modulus increases as the cube of A as a result of 
the mean- field treatment, the zero-shear viscosity increases as A 4 near the gelation point. Obtained dynamic shear 
moduli as a function of the polymer concentration were found to agree well with the experimental data observed for 
aqueous solutions of telechelic poly (ethylene oxide). 

We assumed in this theoretical model that intramolecular associations generating looped chains are absent. Looped 
chains are supposed to compete with intcrmolccular association that causes bridge chains at a junction due to the 
limitation of the multiplicity that the junction can take. Such competition might influence the viscoelasticity of the 
system. This effect as well as the influence of additives such as surfactant or single end-capped polymers will be 
studied in the forthcoming paper. 
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